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Abstract 

We describe a search for evidence of minimal supergravity (mSUGRA) in 92.7 pb _1 of data 
collected with the D0 detector at the Fermilab Tevatron pp collider at ^J~s~ — 1.8 TeV. Events with 
a single electron, four or more jets, and large missing transverse energy were used in this search. The 
major backgrounds are from VF+jets, misidentified multijet, it, and WW production. We observe 
no excess above the expected number of background events in our data. A new limit in terms of 
mSUGRA model parameters is obtained. 



I. INTRODUCTION 

The standard model (SM) has been a great achievement in particle physics. A large number of experimental 
results have confirmed many features of the theory to a high degree of precision. However, the SM is theoretically 
unsatisfactory, and it poses many questions and problems The most notable ones are the fine-tuning problem 

of the SM Higgs self-interaction through fermion loops (^] and the unknown origin of electroweak symmetry breaking 
(EWSB). Supersymmetry (SUSY) B incorporates an additional symmetry between fermions and bosons, and offers 
a solution to the fine-tuning problem and a possible mechanism for EWSB. 

SUSY postulates that for each SM degree of freedom, there is a corresponding SUSY degree of freedom. This results 
in a large number of required supersymmetric particles (sparticles) , and at least two Higgs doublets in the theory. A 
new quantum number, called i?-parity la], is used to distinguish between SM particles and sparticles. All SM particles 
have i?-parity +1 and sparticles have i?-parity —1. The simplest extension to the SM, the minimal supersymmetric 
standard model (MSSM), respects the same SU(3) (£> SU(2) <g) U(l) gauge symmetries as does the SM. SUSY must 
be a broken symmetry. Otherwise we would have discovered supersymmetric particles of the same masses as their 
SM partners. A variety of models have been proposed for SUSY breaking. One of these, the minimal supergravity 
(mSUGRA) model, postulates that gravity is the communicating force from the SUSY breaking origin at a high mass 
scale to the electroweak scale, which is accessible to current high energy colliders. The mSUGRA model is described 
in detail in Ref. g], ft can be characterized by four parameters and a sign: a common scalar mass (mo), a common 
gaugino mass (TO1/2), a common trilinear coupling value (Aq), the ratio of the vacuum expectation values of the two 
Higgs doublets (tan/3), and the sign of //, where n is the Higgsino mass parameter. 

In this analysis, i?-parity is assumed to be conserved. This implies that sparticles must be pair-produced in pp 
collisions. The sparticles can decay directly, or via lighter sparticles, into final states that contain SM particles and 
the lightest supersymmetric particles (LSPs), which must be stable. Because the LSP interacts extremely weakly, it 
escapes detection and leaves a large imbalance in transverse energy ($ T ) in the event. We assume that the lightest 
neutralino (xl) is the LSP, and that Aq = and fi < 0. We fix tan/? = 3 and perform the search in the m 1( /2-mo 
plane. 

Most recently, searches for mSUGRA signatures have been performed at LEP and the Tevatron. At D0, 
dilepton+^ T [Q and jets+$ T || final states have been examined for possible mSUGRA effects. This report de- 
scribes a search in the final state containing a single isolated electron, four or more jets, and large j£ T . One of the 
possible mSUGRA particle-production processes which results in such a final state is shown in Fig. [y. The search is 
particularly sensitive to the moderate too region where charginos and neutralinos decay mostly into SM W and/or Z 
bosons which have large branching fractions to jets. It also complements our two previous searches since the signatures 
are orthogonal to one another. 

II. THE D0 DETECTOR 

DO is a multipurpose detector designed to study pp collisions at the Fermilab Tevatron Collider. The work presented 
here is based on approximately 92.7 pb _1 of data recorded during the 1994-1996 collider runs. A full description of 
the detector can be found in Ref. Here, we describe briefly the properties of the detector that are relevant for this 
analysis. 



3 




FIG. 1. Feynman diagram for gluino pair production and decay to an electron, multijets, and produce $ T . The three-body 
decays are in fact cascade decays in which off-shell particles or sparticles are produced. 



The detector was designed to have good electron and muon identification capabilities and to measure jets and 
ifjrp with good resolution. The detector consists of three major systems: a non-magnetic central tracking system, a 
uranium/liquid-argon calorimeter, and a muon spectrometer. A cut-away view of the detector is shown in Fig. ||. 

The central detector (CD) consists of four tracking subsystems: a vertex drift chamber, a transition radiation 
detector, a central drift chamber, and two forward drift chambers. It measures the trajectories of charged particles 
and can discriminate between singly-charged particles and e + e~ pairs from photon conversions through the ionization 
measured along their tracks. It covers the pseudorapidity region | r7<i j < 3.2. 

The calorimeter is divided into three parts: the central calorimeter (CC) and the two end calorimeters (EC), each 
housed in its own steel cryostat, which together cover the pseudorapidity range |77d| < 4.2. Each calorimeter consists 
of an inner electromagnetic (EM) section, a fine hadronic (FH) section, and a coarse hadronic (CH) section. Between 
the CC and the EC is the inter-cryostat detector (ICD), which consists of scintillator tiles. The EM portion of the 
calorimeters is 21 radiation lengths deep and is divided into four longitudinal segments (layers). The hadronic portions 
are 7-9 nuclear interaction lengths deep and are divided into four (CC) or five (EC) layers. The calorimeters are 
segmented transversely into pseudoprojective towers of Arj x A<fi — 0.1 x 0.1. The third layer of the EM calorimeter, 
where most of the EM shower energy is expected, is segmented twice as finely in both r] and <f>, with cells of size 
A77 x Acj) = 0.05 x 0.05. The energy resolution for electrons is a{E)/E = 15%/y / E(GeV) © 0.4%. For charged pions, 
the resolution is 50%/y/E(GeV) and for jets 80%/ y/E(GeV). The resolution in $ T is 1.08 GeV + 0.019 ■ X).E T (GeV), 
where Et is the scalar sum of the transverse energies in all calorimeter cells. 

The wide angle muon system (WAMUS), which covers \rj^\ < 2.5, is also used in this analysis. The system consists 
of four planes of proportional drift tubes in front of magnetized iron toroids with a magnetic field of 1.9 T and two 
groups of three planes of proportional drift tubes behind the toroids. The magnetic field lines and the wires in the drift 
tubes are transverse to the beam direction. The muon momentum p is measured from the muon's angular bend in the 
magnetic field of the iron toroids, with a resolution of er(l/p) = 0.18(p - 2 GeV)/p 2 © 0.003 GeV -1 , for p > 4.0 GeV. 

A separate synchrotron, the Main Ring, lies above the Tevatron and goes through the CH calorimeter. During 
data-taking, it is used to accelerate protons for antiproton production. Particles lost from the Main Ring can deposit 
significant energy in the calorimeters, increasing the instrumental background. We reject much of this background at 
the trigger level by not accepting events during beam injection into the Main Ring, when losses are largest. 
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Muon Chambers 




FIG. 2. Cut away isometric view of the D0 detector. 
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III. EVENT SELECTION 



Event selection at D0 is performed at two levels: online selection at the trigger level and offline selection at 
the analysis level. The algorithms to reconstruct the physical objects (electron, muon, jet, $ T ) as well as their 
identification at the online and offline levels are described in Ref . Jll| . We summarize below the selections pertaining 
to this analysis. 

A. Triggers 

The D0 trigger system reduces the event rate from the beam crossing rate of 286 kHz to approximately 3-4 Hz, at 
which the events are recorded on tape. For most triggers (and those we use in this analysis) we require a coincidence 
in hits between the two sets of scintillation counters located in front of each EC (level 0). The next stage of the 
trigger (level 1) forms fast analog sums of the transverse energies in calorimeter trigger towers. These towers have a 
size of A?/ x Ac/) = 0.2 x 0.2, and are segmented longitudinally into EM and FH sections. The level 1 trigger operates 
on these sums along with patterns of hits in the muon spectrometer. A trigger decision can be made between beam 
crossings (unless a level 1.5 decision is required, as described below). After level 1 accepts an event, the complete 
event is digitized and sent to the level 2 trigger, which consists of a farm of 48 general-purpose processors. Software 
filters running in these processors make the final trigger decision. 

The triggers are defined in terms of combinations of specific objects required in the level 1 and level 2 triggers. 
These elements are summarized below. For more information, see Refs. P,pT[. 

To trigger on electrons, level 1 requires that the transverse energy in the EM section of a trigger tower be above a 
programmed threshold. The level 2 electron algorithm examines the regions around the level 1 towers that are above 
threshold, and uses the full segmentation of the EM calorimeter to identify showers with shapes consistent with those 
of electrons. The level 2 algorithm can also apply an isolation requirement or demand that there be an associated 
track in the central detector. 

For the later portion of the run, a "level 1.5" processor was also available for electron triggering. In this processor, 
each EM trigger tower above the level 1 threshold is combined with the neighboring tower of the highest energy. The 
hadronic portions of these two towers are also combined, and the ratio of EM transverse energy to total transverse 
energy in the two towers is required to be > 0.85. The use of a level 1.5 electron trigger is indicated in the tables 
below as an "EX" tower. 

The level 1 muon trigger uses the pattern of drift tube hits to provide the number of muon candidates in different 
regions of the muon spectrometer. A level 1.5 processor can also be used to put a px requirement on the candidates 
(at the expense of slightly increased dead time). At level 2, the fully digitized event is available, and the first stage 
of the full event reconstruction is performed. The level 2 muon algorithm can also require the presence of energy 
deposition in the calorimeter consistent with that from a muon. 

For a jet trigger, level 1 requires that the sum of the transverse energies in the EM and hadronic sections of a 
trigger tower be above a programmed threshold. Level 2 then sums calorimeter cells around the identified towers (or 
around the -E/r-weighted centroids of the large tiles) in cones of a specified radius AR = y/ Arf + A(f> 2 , and imposes 
a threshold on the total transverse energy. 

The Iprp in the calorimeter is computed both at level 1 and level 2. For level 1, the vertex z position is assumed to 
be at the center of the detector, while for level 2, the vertex z position is determined from the relative timing of hits 
in the level scintillation counters. 

The trigger requirements used for this analysis are summarized in Table |. Runs taken during 1994-1995 
(Run lb) and during the winter of 1995-1996 (Run lc) were used, and only the triggers "ELE_JET_HIGH" and 
"ELE.JET.HIGHA" in the table were used to conduct this search for mSUGRA. The "EM1.EISTRKCC.MS" trigger 
was used for background estimation. As mentioned above, these triggers do not accept events during beam injection 
into the main ring. In addition, we do not use events which were collected when a Main Ring bunch passed through the 
detector or when losses were registered in monitors around the Main Ring. Several bad runs resulting from hardware 
failure were also rejected. The "exposure" column in Table | takes these factors into account. 

B. Object Identification 
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TABLE I. Triggers used during Run lb and Run lc. 
taking into account the Main Ring vetoes and bad runs. 



'Exposure" gives the effective integrated luminosity for each trigger, 



Trigger Name 



EM1.EISTRKCC.MS 



Exposure 
(Pb- 1 ) 



Level 1 



Level 2 



Run 
period 



82.9 



1 EM tower, E T > 10 GeV 
1 EX tower, E T > 15 GeV 



1 EM tower, E T > 12 GeV, \rj\ < 2.6 

2 jet towers, E T > 5 GeV, r; < 2.0 



1 isolated e, E T > 2Q,GcV 
ffir al > 15 GeV 



Run lb 



ELE..IET.HIGH 



82.9 



1 e, E T > 15 GeV, |r)T< 2.5 
2 jets (Afl, = 0.3), E T > 10 GeV, |r)| < 2.5 
> 14 GeV 



Run lb 



ELE..IET.HIGH 



0.89 



ditto 



ditto 



1 EM tower, E T > 12 GeV, \rj\ < 2.6 

2 jet towers, E T > 5 GeV, |r;| < 2.0 



Run lc 



ELE..IET.HIGHA 



8.92 



1 e, E T > 17 GeV, |r)| < 2.5 
2 jets {AR = 0.3), E T > 10 GeV, \r)\ < 2.5 
#^ al > 14 GeV 



Run lc 



a ^j? is the missing Et in the calorimeter, obtained from the sum of transverse energy of all calorimeter cells. lp T is the 
missing Et corrected for muon momentum, obtained by subtracting the transverse momenta of identified muons from ^? . 



1. Electrons 



Electron identification is based on a likelihood technique. Candidates are first identified by finding isolated clusters 
of energy in the EM calorimeter with a matching track in the central detector. We then cut on a likelihood constructed 
from the following five variables: 

• a x 2 from a covariance matrix that checks the consistency of the shape of a calorimeter cluster with that expected 
of an electron shower; 

• an electromagnetic energy fraction, defined as the ratio of the portion of the energy of the cluster found in the 
EM calorimeter to its total energy; 

• a measure of consistency between the trajectory in the tracking chambers and the centroid of energy cluster 
(track match significance); 

• the ionization deposited along the track dE/dx] 

• a measure of the radiation pattern observed in the transition radiation detector (TRD). (This variable is used 
only for CC EM clusters because the TRD does not cover the forward region ||.) 

To a good approximation, these five variables are independent of each other. 

High energy electrons in mSUGRA events tend to be isolated. Thus, we use the additional restriction: 

£em(0.2) < °' ' (3 ' 



where E tot (0.4) is the energy within AR < 0.4 of the cluster centroid (AR — \J Arf + A</> 2 ) and £'em(0.2) is the 
energy in the EM calorimeter within AR < 0.2. We denote this restriction the "isolation requirement." 

The electron identification efficiency, ef d , is measured using the Z — > ee data. Since only CC (\r) d \ < 1.1) and 
EC (1.5 < < 2.5) regions are covered by EM modules, electron candidates are selected and their identification 
efficiencies are measured in these two regions. An electron is considered a "probe" electron if the other electron in the 
event passes a strict likelihood requirement. This gives a clean and unbiased sample of electrons. We construct the 
invariant mass spectrum of the two electron candidates and calculate the number of background events, which mostly 
come from Drell-Yan production and misidentified jets, inside a Z boson mass window. After background subtraction, 
the ratio of the number of events inside the Z boson mass window before and after applying the likelihood and isolation 
requirements to each probe electron, gives ef d . 

The e? d is a function of jet multiplicity in the event. The presence of jets reduces ef d primarily due to the isolation 
requirement and reduced tracking efficiency. However, with a larger numbers of jets (> 3) in the event, the efficiency 
of locating the correct hard-scattering vertex increases. The two effects compensate each other for events with high 
jet multiplicity [Q . The electron identification efficiencies used in this analysis are obtained from Z — > ee data with 
at least two jets and are given in Table || 

Sometimes a jet with very similar characteristics to an electron can pas s the electron identification selection, and 



result in a fake electron. The effect of fake electrons is discussed in section VA 
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TABLE II. Electron ID efficiencies used in this anaiysis. 



Detector Region CC EC 



0.674 ± 0.039 0.242 ± 0.075 



TABLE III. Parameters for jet identification efficiency as defined in Eq. 



Fiducial Region 


E J T (GcV) 


Po 


pi (GcV" 1 ) 


p 2 (GcV" 2 ) 


CC 


15-27.4 


0.8994 ± 0.0070 


(5.04 ± 0.45) x 10 - ' 1 


(-6.7 ± 1.0) x 10- a 


< i-o) 


> 27.4 


0.9864 ± 0.0005 


(2.16 ± 0.57) x 10 -5 


(-1.90 ± 0.30) x 10~ 7 


ICR 


15-30.5 


0.9838 ± 0.0017 


(9.76 ± 1.33) x 10 -4 


(-1.76 ± 0.27) x 10~° 


(i.o <\ni\< is) 


> 30.5 


0.9981 ± 0.0008 


(-2.27 ± 2.26) x 10~ 5 


(-1.52 ± 1.22) x 10~ 7 


EC 




0.9866 ± 0.0004 


(-3.81 ± 1.05) x 10~ 5 


(-1.15 ± 0.75) x 10~ y 


(1.5 < \rf d \ < 2.5) 











2. Jets 



Jets are reconstructed in the calorimeter using a fixed-size cone algorithm with AR = 0.5. A jet that originates 
from a quark or a gluon deposits a large fraction of its energy in the FH part of the calorimeter, and so we identify 
jets through the fractional energy in the EM and CH parts of the calorimeter. We require the fraction of the total 
jet energy deposited in the EM section of the calorimeter (em/) to be between 0.05 and 0.95 for high energy jets 
(E? T > 35 GeV), and the fraction of the total jet energy deposited in the CH section of the calorimeter (chf ) to be less 
than 0.4. Because electronic and uranium noise is generally of low energy, the lower bound of the emf requirement 
is raised gradually for lower energy jets in the CC. (It is 0.2 for CC jets with E 3 T k, 15 GeV.) Because there is no 
electromagnetic coverage in the ICR, we do not apply a lower bound cut on emf in that region. 

A multijet data sample corrected for detector noise is used to measure the jet identification efficiency, e\ d . The 
efficie ncy is a function of E 3 T , and is parametrized as in Eq. 3.2, with the fitted values of the parameters listed in 



Table III. 



i d =Po+PixE J T +p 2 x {E 3 T f. (3.2) 



3. Muons 

To avoid overlapping with the dilepton analysis, we veto events containing isolated muons satisfying all the following 
criteria: 

• The muon has a good track originating from the interaction vertex. 

• The muon has pseudorapidity \rj^. \ < 2.5. 

• There is a large integrated magnetic field along the muon trajectory (J B ■ dl ). This ensures that the muon 
traverses enough of the field to give a good Pt measurement. 

• The energy deposited in the calorimeter along a muon track is at least that expected from a minimum ionizing 
particle. 

• Transverse momentum pt > 4 GeV. 

• The distance in the r\ — <f> plane between the muon and the closest jet is AR(p,j) > 0.5. 

4- Event selection 

About 1.9 million events passed the ELE_JET_HIGH and the ELE_JET_HIGHA triggers. We require at least 
one electromagnetic cluster with Et > 18 GeV and a track matched to it. The interaction vertex must be within 
\z v \ < 60 cm. About 600,000 events remain after these selections. Kinematic and fiducial requirements are then 
applied to select our base data sample. The criteria are listed below, with numbers in the curly brackets specifying 
the number of events surviving the corresponding requirement. 
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• One electron in the good fiducial volume (|?7|| < 1.1 or 1.5 < 1 77^ | < 2.5) passing restrictive electron identification 
criteria, and with EJ, > 20 GeV — {15547}. 

• No extra electrons in the good fiducial volume passing "loose" electron identification for > 15 GeV. The 
selection criteria for the "loose" electrons are the same as those used for signal electrons in the dilepton analysis, 
keeping two analyses independent of each other — {15319}. 

• |?f I < 2.0 — {13997}. 

• No isolated muons — {13980}. 

• Four or more jets with E ] T > 15 GeV and \rf A \ < 2.5 — {187}. 

• $ T > 25 GeV — {72}. 

After these selections the base sample contains 72 events. The major SM backgrounds are from W+ > 4 jets 
— > e + v+ > 4 jets, tt — > Wb Wb — > e + v+ > 4 jets, WW+ > 2 jets — > e + v+ > 4 jets, and multijet events in which 
one of the jets is misidentified as an electron and the jet transverse energies are inaccurately measured to give rise to 



IV. EVENT SIMULATION 

We use PYTHIA |l3| to simulate mSUGRA signal and tt and WW backgrounds. We check our results and obtain 
generator-dependent systematic errors using the herwig Q generator. W boson and associated jet production is 
generated using VECBOS jl5| and herwig. The final state partons, which are generated by vecbos as a result of a 
leading order calculation, are passed through herwig to include the effects of additional radiation and the underlying 
processes, and to model the hadronization of the final state partons Jig] . 

In order to efficiently search for mSUGRA in a large parameter space and to reduce the statistical error on signal 
acceptance, we used a fast Monte Carlo program called FMC0 to model events in the D0 detector and to calculate 
the acceptance for any physics process passing our trigger and offline selections. The flow-chart of FMC0 is shown 
in Fig. ||. First, through a jet-reconstruction program, the stable particles that interact in the detector are clustered 
into particle jets, in a way similar to the clustering of calorimeter cells into jets. However, the generated electrons, if 
they are not close to a jet (AR > 0.5 in 77 — <j> space), are considered as the electrons reconstructed in the detector. 
Otherwise, they are clustered into the jet. The generated muons are considered as the reconstructed muons in the 
detector. Next, the electrons, jets, muons, a nd $ T in the events are smeared according to their resolutions determined 



from data [|l6|. The offline selections (Sec. [II B 4 ) are applied to the smeared objects. Finally, each passed event is 
weighted with trigger and identification efficiencies. The outputs of FMC0 are an "ntuple" that contains the kinematic 
characteristics (Et, f], 4>, etc.) of every object and a run-summary ntuple that contains the information of trigger 
efficiency and total acceptance for the process being simulated. The acceptance A is calculated as follows: 

N p „ 

^ = ^-E^'-^-4 ts , (4.1) 

JV gen . 

where sl°^ 1 is the overall trigger efficiency, ef d is the electron identification efficiency, e^ s is the product of jet 
identification efficiencies of the four leading jets, N gen is the number of generated events, and iV pass is the number of 
events that pass the offline kinematic requirements. The uncertainty on the acceptance, Sa, is calculated as: 

iVgen . 

where 5 £ comes from the propagation of uncertainties on ej^g 1 , ef d , and £ J id ts . Since the same electron and jet 
identification efficiencies, and t he same trigger turn-ons are used the error on the acceptance is 100% correlated 



event-by-event as shown in Eq. 4.2 



Because the signal triggers impose a combination of requirements on the electron, jets, and $ T , the overall trigger 
efficiency has three corresponding components. The efficiency of each component was measured using data. The 
individual efficiencies are then used to construct the overall trigger efficiency. The details of the measurements and 



construction are documented in Ref. 12 1. Table |TV| compares the trigger efficiencies of W+ jets events measured 
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FIG. 3. Flow-chart of FMC0. Prefix "s_" refers to smeared objects. 



TABLE IV. Comparison of e^lf, the total trigger efficiency of ELE_JET_HIGH trigger. The second column lists the 
efficiencies measured using W+ jets data; the third column lists the simulated efficiencies found by putting the VECBOS 
W+ jets events through FMC0. 





Data 


VECBOS 


> 1 


0.589 ±0.019 


0.579 ± 0.022 


> 2 


0.826 ± 0.027 


0.833 ± 0.020 


> 3 


0.928 ±0.031 


0.925 ±0.016 


> 4 


0.944 ± 0.037 


0.957 ±0.012 









in data with those simulated using VECBOS Monte Carlo. We find that they are in good agreement at each jet 
multiplicity. 

We also compared the acceptance of FMC0 with GEANT Jl8| and data, and found good agreement for W+ jets, ti, 
and WW events. 



V. BACKGROUNDS 



A. Multijet background 

From the ELE_JET_HIGH and ELE_JET_HIGHA triggered data we obtain two sub-samples. For sample 1, we 
require all offline criteria to be satisfied, except for $J T . At small $ T (< 20 GeV), sample 1 contains contributions 
mainly from multijet production, where jet energy fluctuations give rise to $ T . At large $ T (> 25 GeV), it has 
significant contributions from W+ jets events, with additional contributions from ti production and possibly the 
mSUGRA signal. For sample 2, we require that the EM object represent a very unlikely electron candidate by 
applying an "anti-electron" requirement |l2| |. All other event characteristics are the same as those in sample 1. The 
sample 2 requirements tend to select events in which a jet mimics an electron, and consequently sample 2 contains 
mainly multijet events with little contribution from other sources for Jj] T > 25 GeV. The $ T spectra of the two samples 
can therefore be used to estimate the number of multijet background events (A^muitijet) in sample 1 as follows. We 
first normalize the $J T spectrum of sample 2 to that of sample 1 in the low-$ T region, and then estimate A^uitijet by 
multiplying the number of events in the signal region (Jf, T > 25 GeV) of sample 2 by the same relative normalization 
factor [19]. 

The lf]rp spectra for both samples are shown in Fig. ||, normalized to each other for < ]p T < 14 GeV, and for the 
cases in which the fake electron is in the CC and EC, respectively. From these distributions, we calculate iV mu itijot to 
be 82.6 ± 15.3 and 19.1 ± 4.7, for inclusive jet multiplicities of 3 and 4 jets, respectively. (The inclusive 3-jet sample 
is obtained the same way as the base sample, except that we require at least 3 jets, rather than 4, in the event.) The 
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errors include statistical uncertainties and systematic uncertainties in the trigger and object identification efficiencies, 
different definitions of sample 2, and different choice for the normalization regions. 



B. tt background 

The number of background events, N t j, is calculated using FMC0. The it events were generated using PYTHIA E3] 
for m top = 175 GeV. A tt production cross section of a = 5.9 ± 1.7 pb, as measured by D0 |pof , is used. The results 
are N t j = 27.7 ± 8.3 events and N t j = 16.8 ±5.2 events for inclusive jet multiplicities of 3 and 4 jets, respectively. The 
errors include uncertainties on the tt production cross section, differences in physics generators, trigger and object 
identification efficiencies, and on the integrated luminosity. 



C. WW+ jets background 



FMC0 is also used to calculate the WW+ jets background. The production cross section at next-to-leading order 
is taken as a = 10.40 ± 0.23 pb J2l],g2|], assuming no anomalous couplings (5n = A = 0) |23|j . The WW events were 
generated using PYTHIA. There are 7.7 ± 1.2 and 1.4 ± 0.3 events expected for inclusive jet multiplicity of 3 and 4 
jets, respectively. The errors include uncertainties on the production cross section, trigger and object identification 
efficiencies, differences in physics generators, the jet energy scale, and on the integrated luminosity. 



D. W+ jets background 

To good approximation, each extra jet in W+ jets events is the result of an extra coupling of strength a s and 
we expect the number of W+ jets events to scale as a power of Aj et . The scaling law is supported by the W+ jets, 
Z+ jets, and 7+ jets data ^4|. In this analysis, we first estimate the number of W+ > 3-jet events, , in the data 
collected with ELE_JET_HIGH and ELE_JET_HIGHA triggers, and then extract the effective scaling factor a using 
W+ > n-jet events collected with EM1_EISTRKCCJVIS trigger. The expected number of W+ > 4-jet events (N^) 
in our base sample is then: 

Wi 

N w = N w. a .«i t (5 . 1} 

^trig 



where gffi and are trigger efficiencies of W+ > 3-jet and W+ > 4-jet events, respectively, as shown in Table IV 



1. Estimating the number ofW+ > n-jet events 

We estimate the number of W+ > n-jet events in a way similar to that used to estimate the multijet background. We 
first use a neural network (NN) to define a kinematic region in which W+ > n-jet events dominate the background and 
any possible contribution from mSUGRA can be neglected. In that region, we normalize the number of W+ > n-jet 
MC events to the number of events observed in the data which have had all other major SM backgrounds subtracted. 
The normalization factor is then applied to the whole W+ > n-jet MC sample to obtain our estimate for the W+ > n- 
jet background in the data. 

In this analysis, we use a NN package called MLPFIT [E5[. All NNs have the structure of X-2X-1, where X is the 
number of input nodes, i.e., the number of variables used for training, and 2X is the number of nodes in the hidden 
layer. We always use 1 output node with an output range of to 1. Signal events (in this case, W+ > n-jet events) 
are expected to have NN output near 1 and background events near 0. We choose the NN output region of 0.5-1.0 
to be the "signal" -dominant kinematic region. The variables used to distinguish W+ > n-jet events from other SM 
backgrounds and the mSUGRA signal are: 

• $T 

• E% 

• H T = J2 E T for a11 J etS with E T > 15 GeV 
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CC (N= t >4) 




100 



20 40 60 80 100 



FIG. 4. $ T spectra of sample 2 (points) normalized to sample 1 (histograms) in the region of < $ T < 14 GeV. The 
normalizations are done for the fake electron in the CC and EC, respectively. The errors are statistical only. 
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• &4>e,$ T 

• M T = ^2^^ T [l-cos(A0 e ^ r )] 

• A(f)j 1> f< T (not used for > 4-jet events) 

• A(j>j 2 ^ (used for > 2-jet and > 3-jet events) 

• A — aplanarity pr| (used for > 2, > 3, and > 4-jet events) is denned in terms of the normalized momentum 
tensor of the W boson and the jets with E 3 T > 15 GeV: 

M ab = E if g f , (5.2) 

where pi is the three-momentum of object i in the laboratory frame, and a and b run over the x, y, and z 
coordinates. Denoting Qi, Q2, and Q3 as the three eigenvalues of M t in ascending order, A — 1.5 x Q\. The 
p z of the W boson is calculated by imposing the requirement that the invariant mass of the electron and the 
neutrino (assumed to be the source of f T ) equals the W boson mass. This requirement results in a quadratic 
equation for the longitudinal momentum of the neutrino. Because the probability of a small p z is usually higher 
than that of a large p Zl the smaller p z solution is always chosen. In cases where there is no real solution, f T is 
increased until a real solution is obtained. 

• r# = Ht2/Hz, where H?2 — Ht — Ei£ , and Hz — £V \Pz\ wnere i runs over the electron, all jets with 
Eip > 15 GeV, and neutrino (as assumed in the calculation of A) in the event (27) (only used for > 4-jet events). 

• cos 9*, where 9* is the polar angle of the electron in the W boson rest frame, relative to the direction of flight 
of the W boson. The W boson four-momentum is obtained by fitting the event to a ti assumption. The details 
of the fit are described in Ref. J2?j (only used for > 4-jet events). 

• cos 9* b , where 8* eb is the angle between the electron and the b jet from the same top (or antitop) quark in the W 
boson rest frame p8| . Again, a fit to the tt assumption is performed to identify the correct b jet (only used for 
> 4-jet events). 



All the offline requirements described in Sec. [II B 4 are applied except that the requirement on the number of 



jets is reduced corresponding to different inclusive jet multiplicity. The multijet, tt, and WW backgrounds are 



estimated using the methods described in Sec. VA-VD. The mSUGRA events were generated with m — 170 GeV, 
m i/2 — 58 GeV and tan/3 = 3. This parameter set was chosen because it is close to the search limit obtained in the 
dilcpton analysis. 



2. Estimating N3 

The result of the NN training for > 3-jet events is shown in Fig. |^(a). The number of W+ > 3-jet events used in the 
training is the same as the sum of all background events, including any possible mSUGRA sources in their expected 
proportions. The match between training and data is shown in Fig. ^(b), where the data and MC are normalized to 
each other for NN output between 0.5 and 1.0. Because the number of mSUGRA events is negligible in this region, 
we do not include them in the background subtraction. We estimate that 241.8 ± 18.0 W+ > 3-jet events pass our 
final 3-jet selection. 



3. Measuring the scaling factor a 



We extract the parameter a from the data passing the EM1_EISTRKCC_MS trigger, which does not have a jet 

w 

requirement in the trigger, and fit the measured number of W + n-jet events (N n ) to: 



— w -^w 



(5.3) 



7"" 



N n values are obtained as described in Sec. VD 1. The NN training and normalization to the data are performed 

. . r— 1 — w 

separately for each inclusive jet multiplicity. The results are summarized in Table M. The errors on N n include 
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(b) 



FP 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



NN Output 
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NN Output 



FIG. 5. (a) Expected NN output for events passing the ELE_JET_HIGH or ELE_JET_HIGHA triggers and with > 3 jets, 
(b) Expected NN output for data (points) and the observed NN output for data (histogram). The error on the points include 
statistical and systematic errors. All events were required to pass our offline selections, except that we required only 3 jets 
instead of 4. 

— w 

TABLE V. Estimated number of W+ > n-jet events, N n , as a function of inclusive jet multiplicity in the data passing the 
EM1_EISTRKCCJVIS trigger. They were obtained by normalizing MC to data in the NN output region where W+ > n-jets 
events dominate (see text). iVdata is the number of observed events. The mSUGRA events were generated with mo = 170 GeV, 
m i/2 = 58 GeV, and tan f3 — 3. 



> 1 



> 2 



> 3 



> 4 



■IVdata 
-^multijet 

N t -t 

Nww 
— w 
N 

1 v n 

NmSVGRA 



8191 
826 ± 95 
25.8 ± 7.6 
33.7 ±3.3 
7210 ± 131 
28.3 ±3.7 



1691 
291 ± 48 
26.1 ± 7.6 
23.6 ±2.3 
1283 ± 79 
25.0 ±3.1 



353 
75 ± 15 
21.9 ±6.5 
6.19 ±0.95 
230 ± 27 
19.7 ± 2.7 



64 

16.6 ± 7.0 

13.5 ±4.3 
1.12 ±0.25 
27.4 ± 7.4 

12.6 ±2.1 



statistical errors from MC and data, and uncertainties on the choice of different normalization regions and on the 
choice of different QCD dynamic scales used in generating vecbos events. 

— W i 1 in 

The fit of N n to Eq. |5J is shown in Fig. § from which we extract a = 0.172 ± 0.007. 



4- Calculating the number ofW+ > 4-jef events, N4 



With = 0.925 ± 0.016 and eg£ = 0.957 ± 0.012, and using Eq. 5.1, we obtain N? = 43.0 ± 7.6. 



E. Summary 



The expected numbers of events in the base data sample from the major sources of background are summarized in 
Table VI. From the table, we conclude that the sum of the backgrounds is consistent with the observed number of 
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NY = NY • a 

a = 0.172 ±0.007 
% 2 /dof = 1.1 



w 



i-1 



3 4 5 

Inclusive jet multiplicity 



FIG. 6. Fit of W+ > n-jet events to the power law of Eq. 5.3 
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TABLE VI. Expected numbers of events in the base data sample 
observed data events. 


from the major sources of background and the number of 


W+ > 4-jets 


43.0 ± 7.6 


misidentified multijet 


19.1 ± 4.7 


it 


16.8 ± 5.2 


WW+ > 2-jets 


1.4 ± 0.3 


Total 


80.3 ± 10.4 


Data 


72 



candidate events. 



VI. SEARCH FOR SIGNAL 
A. Neural Network Analysis 



We use a NN analysis to define a kinematic region in which the sensitivity of signal t o back ground is highest. We 



use the following variables in the NN. Those not defined below have been defined in Sec. |VD 1 



• ]fjrp - For the signal, $ T comes from two LSPs and at least one neutrino. For the it, W+ jets, and WW 
backgrounds, it comes from the neutrino. For multijet background, it comes from fluctuation in the measurement 
of the jet energy. Generally, the signal has larger $ T than the backgrounds. 

• E T - The electron in the signal comes from a virtual W boson decay. Its spectrum is softer than that of the 
electrons from the it and W+ jets backgrounds. 

• Ht - A pair of heavy mSUGRA particles are produced in the hard scattering and most of the transverse energy 
is carried away by jets. The Ht for the signal thus tends to be larger than that for the major backgrounds. 

• Ej? - The third leading jet in Et from W+ jets, WW, and multijet events most likely originates from gluon 
emission. For it and mSUGRA events, it is probably due to W boson decay. Thus, the it and mSUGRA signals 
have a harder Eij? spectrum. 

• Mt - For it, W+ jet, and WW events, Mt peaks near Mw = 80 GeV. This is not the case for the signal since 
we expect the W boson produced in the decay chain to be virtual for a wide range of mx/i up to 200 GeV. 

• Acj) e j£ T - Because the electron and neutrino form a W boson in it, W+ jet, and WW events, their Ac/> e ^ 
spectra should peak away from A0 e ^ = 0. For multijet events, the A(j> e ^ spectrum should peak near and 
7r because ~$ T can be caused by fluctuations in the energy of the jet which mimics an electron. 

• A - W+ jets, WW, and multijet events are more likely to be collinear due to QCD bremsstrahlung, while the 
signal and it events are more likely to be spherical. 

• cos 9*, where 9* is the polar angle of the higher-energy jet from W boson decay in the rest frame of parent W 
boson, relative to the direction of flight of the W boson. This is calculated by fitting all the events to the it 
assumption. For it production, the spectrum is isotropic, but for the signal and other SM backgrounds, it is 
not. 

• cos 9*, the signal has a somewhat different cos#* distribution than the background does, especially for it events. 

The spectra for these variables are shown in Fig. [?]. There is no evidence of an excess in our data for the mSUGRA 
parameters used. Fig. || displays the cos 9* and cos 9* distributions for signal and it events. These two variables are 
particularly useful in reducing the it background relative to the mSUGRA signal. Nevertheless, it events still make 
the largest contribution in the signal-rich region because of their similarity to the mSUGRA signal. This can be seen 
in Fig. |[ in which the NN output is displayed for each background and the mSUGRA signal for a particular set of 
parameters. The result of the NN output for data is given in Fig. |l0|. The expected background describes the data 
well. 
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FIG. 7. Distribution of NN variables for data (open histogram), background (points) and signal (hatched histogram). The 
signal was generated at mo = 170 GeV, TO1/2 = 58 GeV, and tan/3 = 3. We have multiplied the expected number of signal 
events (18.5) by a factor of 4.3 to normalize it to the total number of background events. Since the same number of signal and 
background events are used to train the NN, the plot shows the relative strength of signal to background as seen by the NN. 
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FIG. 8. Distribution of (a) cos 9* and (b) cos 8* for signal (hatched histogram) and tt events (points). The signal was 
generated at mo = 170 GeV, TO1/2 = 58 GeV, and tan/3 = 3. We have multiplied the expected number of signal events (18.5) 
by a factor of 0.91 to normalize it to the number of tt events expected in our base sample. 



B. Signal Significance 



To apply the optimal cut on the NN output, we calculated the signal significance based on the expected number of 
signal (s) and background (b) events that would survive any NN cutoff. We define the significance (S) below. The 
probability that the number of background events, b, fluctuates to n or more events is: 



F{n\b) 



k—n 



P {k\b) 



i 



'2it 



- t2 / 2 dt 



(6.1) 



S(n\b) 



where p(k\b) = 2-S — is the Poisson probability for observing k events with b events expected. S(n\b) can be regarded 
as the number of standard deviations required for b to fluctuate to n, and it can be calculated numerically. For s + b 
expected events, the number of observed events can be any number between [0, oo). The significance is thus defined 
as: 



S = ^Tp(n\s + b)-S(n\b) 



(6.2) 



n=0 



where p(n\s + b) is the Poisson probability for observing n events with s + b events expected. 

The NN output corresponding to the maximum significance determines our cutoff to calculate the 95% C.L. limit 
on the cross section. The error on the expected signal includes uncertainties on trigger and object identification 
efficie ncies , on parton distribution functions (10%), differences between MCs (12%), and on the jet energy scale (5%). 
Table VII lists the results in terms of 95% C.L. limits on production cross sections for various sets of model parameters 
of mSUGRA. 



VII. RESULTS 



We conduct an independent NN analysis on each generated mSUGRA point. The production cross section calculated 
by PYTHIA is compared with that obtained by limit calculation at 95% C.L. to determine whether the mSUGRA point 
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NN Output 



FIG. 9. Result of training of a NN. The excess above the background near 1 is the expected signal. The signal was generated 
at mo = 170 GeV, mi/2 = 58 GeV, and tan/3 = 3. The backgrounds are stacked up in the order of W(ev)+ jets, W(tv) + jets, 
misidentified multijet, tt, and WW production. The contribution of each type of background is normalized to its expected 
number of events in the data. 



19 




0.8 1 1.2 
NN Output 



FIG. 10. NN output for data (open histogram), signal (hatched histogram), and background (points). The signal was 
generated at mo = 170 GeV, mi/ 2 = 58 GeV, and tan j3 — 3. The background expectation describes the data well. The vertical 
arrow i ndica tes the cutoff on the NN output that corresponds to the maximum signal significance. The significance (described 
in Sec. VI B) as a function of NN output is plotted in the insert. 
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TABLE VII. Number of observed events (N ohB ), expected total background events (AT^gd), and expected signal events 
(A^msuGRA), corresponding to the optimal NN cutoff for different sets of mSUGRA parameters. The signal acceptance after 
NN cutoff (Acceptance), mSUGRA production cross section for each parameter set (er mSUGRA) , and the calculated 95% C.L. 
upper limit on the production cross section (095%) are also listed. All limits are for tan/3 = 3. 



m 


m 1/2 


N ohB 


a rtotal 
Jv bkgd 


A'mSUGRA 


Acceptance 


0" mSUGRA 


°95% 


(GeV) 


(GeV) 








(%) 


(Pb) 


(Pb) 


160 


60 


8 


6.45 ± 1.22 


11.11 ± 1.97 
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is excluded or not. Using the two cross sections at each point, we linearly extrapolate between the excluded and non- 
excluded points to determine the exact location of the exclusion contour. The exclusion contour at the 95% C.L. is 
plotted in Fig. [n|. Shown in the same figure are the results of the D0 dilepton and LEP I |2^] analyses. 

Our single-electron analysis is particularly sensitive in the moderate too region. The extended region of exclusion 
relative to the D0 dilepton result is in the range of 165 GeV < too < 250 GeV. The dominant SUSY process changes 
from gq production at mo = 170 GeV to g pair production at Too = 250 GeV. The limit worsens as too increases 
because the mass difference between \t {x\) an d Xi decreases, resulting in softer electron and jets spectra, and 
consequently reduced acceptance. 

As this work was being completed, a related result Q on searches for mSUGRA in the jets plus missing energy 
channel at Tevatron appeared. Since its limits on mSUGRA parameters, although more restrictive than those obtained 
in this work and in the earlier D0 publication |31j| in the analogous channel, are expressed in a different parameter 
plane (to- vs. to~), we do not show them in Fig. |11|. 



VIII. CONCLUSION 



We observe 72 candidate events for an mSUGRA signal in the final state containing one electron, four or more jets, 
and large $ T in 92.7 pb _1 data. We expect 80.3 ± 10.4 such events from misidentified multijet, ti, W+ jets, and WW 
production. We conclude that there is no evidence for the existence of mSUGRA. We use neural network to select a 
kinematic region where signal to background significance is the largest. The upper limit on the cross section extends 
the previously D0 obtained exclusion region of mSUGRA parameter space. 



ACKNOWLEDGMENTS 



We thank the staffs at Fermilab and collaborating institutions, and acknowledge support from the Department of 
Energy and National Science Foundation (USA), Commissariat a L'Energie Atomique and CNRS/Institut National 
de Physique Nucleaire et de Physique des Particules (France), Ministry for Science and Technology and Ministry for 
Atomic Energy (Russia), CAPES and CNPq (Brazil), Departments of Atomic Energy and Science and Education 
(India), Colciencias (Colombia), CONACyT (Mexico), Ministry of Education and KOSEF (Korea), CONICET and 
UBACyT (Argentina), The Foundation for Fundamental Research on Matter (The Netherlands), PPARC (United 
Kingdom), Ministry of Education (Czech Republic), A. P. Sloan Foundation, NATO, and the Research Corporation. 



* Also at University of Zurich, Zurich, Switzerland. 

' Also at Institute of Nuclear Physics, Krakow, Poland. 

[1] C. Quigg, "Gauge Theories of the Strong, Weak, and Electromagnetic Interactions," Addison- Wesley, Reading, MA, 1983, 
p. 188-189. 

[2] R. Cahn, Rev. Mod. Phys. 68, 951 (1996). 
[3] J. L. Hewett, hep-ph/9810316 (unpublished). 

[4] For a review, see e.g., H. Haber and G. Kane, Phys. Rept. 117, 75 (1985). 
[5] G. R. Farrar and P. Fayet, Phys. Lett. B 76, 575 (1978). 
[6] V. Barger et ah, hep-ph/0003154 (unpublished). 

[7] D0 Collaboration, B. Abbott et ah, Phys. Rev. D Rapid Comm. 63, 091102 (2001). 
[8] D0 Collaboration, B. Abbott et ah, Phys. Rev. Lett. 83, 4937 (1999). 
[9] D0 Collaboration, S. Abachi et al, Nucl. Instr. Methods Phys. Res. A 338, 185 (1994). 
[10] At D0, we define 9 and <f> to be the polar and azimuthal angles of the physical objects, respectively. We define the 
pseudorapidity r\ = — ln(tan#/2). We denote r\& as the pseudorapidity calculated relative to the center of the detector 
rather than relative to the reconstructed interaction vertex. 
[11] D0 Collaboration, S. Abachi et al, Phys. Rev. D 52, 4877 (1995). 
[12] J. Zhou, Ph.D. thesis, Iowa State University, 2001 (unpublished). 

http: / / www-d0.fnal.gov / results / publications_talks /thesis /johnzhou / thesis.ps. 



22 




FIG. 11. Exclusion contour at the 95% C.L. for mSUGRA with tan/3 = 3. The result from the D0 dilepton and LEP I 
analyses are also shown. 



23 



[13] M. Bengtsson and T. Sjostrand, Comp. Phys. Comm. 43, 367 (1987); S. Mrenna, Comp. Phys. Comm. 101, 23 (1997); 

and T. Sjostrand, L. Lonnblad, and S. Mrenna, hep-ph/0108264 (unpublished). 
[14] G. Marchesini et al., Comp. Phys. Comm. 67, 465 (1992). We used v5.9. 
[15] W. Beenakker, F. A. Berends, and T. Sack, Nucl. Phys. B 367, 32 (1991). 
[16] D0 Collaboration, S. Abachi et al., Phys. Rev. Lett. 79, 1197 (1997). 
[17] R. Genik, Ph.D. thesis, Michigan State University, 1998 (unpublished). 

http: / / www-dO.fnal.gov / results / publications_talks /thesis /genik / thesis_lite.ps. 
[18] CERN Program Library Long Writeup W5013, 1994. 

[19] J. Yu, Ph.D. thesis, State University of New York, Stony Brook, 1993 (unpublished). http://www- 

dO.fnal.gov/results/publications_talks/thesis/yu/jae_thesis_final.ps. 
[20] D0 Collaboration, B. Abbott et al, Phys. Rev. D 60, 012001 (1999). 
[21] J. Ohnemus, Phys. Rev. D 44, 1402 (1991). 

[22] D0 Collaboration, B. Abbott et al, Phys. Rev. D Rapid Comm. 63 031101 (2001). 

[23] K. Hagiwara et ah, Nucl. Phys. B 282, 253 (1987). 

[24] D0 Collaboration, S. Abachi et ah, Phys. Rev. Lett. 79, 1203 (1997). 

[25] J. Schwindling and B. Mansoulie, "MLPFIT" , 

http:/ /schwind. home. cern.ch/schwind/MLPfit. html (unpublished). 
[26] V. Barger and R. Phillips, "Collider Physics," Addison- Wesley, Reading, MA, 1987, p. 281. 
[27] D0 Collaboration, B. Abbott et ah, Phys. Rev. D 58, 052001 (1998). 
[28] G. Mahlon and S. Parke, Phys. Lett. B 411, 173 (1997). 

[29] LEPSUSYWG, ALEPH, DELPHI, L3, and OPAL experiments, http://www.cern.ch/LEPSUSY. 
[30] CDF Collaboration, T. Affolder et al, Phys. Rev. Lett. 88, 041801 (2002). 
[31] D0 Collaboration, B. Abbott et al, Phys. Rev. Lett. 83, 4937 (1999). 



24 



